Comment on "Existence of Internal Modes of Sine-Gordon Kinks" 
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In Ref.[l] [Phys. Rev. B. 42, 2290 (1990)] we used a rigorous projection operator collective 
variable formalism for nonlinear Klein-Gordon equations to prove the continuum Sine-Gordon (SG) 
equation has a long lived quasimode whose frequency u)b= 1.004 To is in the continuum just above the 
lower phonon band edge with a lifetime (I/ts ) = 0.0017 Fq. We confirmed the analytic calculations 
by simulations which agreed very closely with the analytic results. In Ref.[3] [Phys. Rev. E. 62, R60 
(2000)] the authors performed two numerical investigations which they asserted "show that neither 
intrinsic internal modes nor quasimodes exist in contrast to previous results." In this paper we 
prove their first numerical investigation could not possibly observe the quasimode in principle and 
their second numerical investigation actually demonstrates the existence of the SG quasimode. Our 
analytic calculations and verifying simulations were performed for a stationary Sine-Gordon soliton 
fixed at the origin. Yet the authors in Ref.[3] state the explanation of our analytic simulations 
and confirming simulations are due to the Doppler shift of the phonons emitted by our stationary 
Sine-Gordon soliton which thus has a zero Doppler shift. 

PACS numbers: 05.45.Yv, 03.50.-z 
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I. INTRODUCTION 

In Ref.[l| we proved the continuum Sine-Gordon (SG) 
equation has a long lived quasimode whose frequency 
from simulation is uig = (1.004 ± 0.001) Fq (where Fq is 
the frequency of the lower band edge in units where the 
speed of sound c = 1) and whose lifetime from simulation 
is (I/ts) = (0.003 ± 0.001) Fo . We used a rigorous pro- 
jection operator collective variable (CV) formalism for 
nonlinear Klein- Gordon equations derived in Ref.J3| to 
calculate the quasimode frequency and lifetime. Our cal- 
culated theoretical values for the frequency and inverse 
lifetime are n = 1.00585 Fq and (l/r^) = 0.0017 Fq which 
agree very well with our simulation values. 

In Ref . [3 the authors performed two numerical investi- 
gations for the SG in which they assert "show that neither 
intrinsic internal modes nor quasimodes exist in contrast 
to previous results" referring to Ref.p^. In Sees. Ill and 
IV we analyze their two numerical investigations in de- 
tail and prove their first numerical investigation could not 
possible observe the quasimode in principle and that their 
second numerical investigations actually observes the SG 
quasimode at the beginning of their simulation. However, 
the length of their system was so short for their long 
observation time, that there were many transversals of 
the system by phonons emitted at different times by the 
soliton, which then reflected from the end of the system 
and then interfered with phonons emitted later. Thus 
each phonon interfered with phonons emitted earlier and 
phonons emitted later which led to a very complicated 
interference pattern. The authors of Ref. 3] concluded 
that the complicated interference pattern was a "proof" 
that SG quasimodes don't exist. However, the correct 
conclusion is that their poorly designed numerical inves- 
tigation was for a time t that was more than ten times 
too long for the length of their system to avoid the irrele- 
vant interferences. In the first 200 seconds just before the 
first emitted phonons reflected off the end of the system 



and returned to the stationary emitting SG soliton, the 
finite lifetime of the quasimode is clearly observable. 

Our analytic calculations and verifying simulations 
were all for a continuum, force free and stationary SG 
soliton i.e., the center of mass of the SG was fixed at 
the origin for all times. However, in their two numerical 
investigations in Ref. j3] the author provided an explana- 
tion of our analytic calculations and verifying simulations 
which was that our phonons were Doppler shifted. Which 
is truly amazing for phonons emitted by a stationary SG 
soliton fixed at the origin. As a result, their two numer- 
ical investigations and their "explanation" of our results 
have absolutely no relevance to the validity of our ana- 
lytic soliton and verifying simulations of our continuum 
stationary and force free SG quasimode. 

In Sec. II we outline the derivation of the exact equa- 
tions of motion for the SG equation. We prove in Sec. Ill 
that the first numerical search for the SG quasimode in 
Ref.JSll could not observe the SG quasimode in principle. 
While in Sec. IV we show that in their second numeri- 
cal investigation the authors of Ref.Q actually observe 
the SG quasimode at the beginning of their simulation. 
However, their simulation was for a time, too large for 
the length of their system. Consequently, they observed 
a complicated interference pattern which was totally ir- 
relevant in the SG quasimode mode that was clearly ob- 
servable at the beginning of their simulation. In Sec. VI 
we present our conclusions and discuss a recent work, 
Ref. 13, which contains a new solution of the SG equa- 
tion by using the inverse transform method and find our 
SG quasimode solution is valid. 



II. EQUATIONS OF MOTION FOR THE SG 
QUASIMODE 

The purpose of this section is to outline the derivation 
of the equations of motion from our exact CV equations 



for the SG equation in Ref.[ll| that we actually solved 
for the SG quasimode in Ref.P]. We need these rigor- 
ous equations in order to contrast them with the two CV 
equations of motion for the kink momentum P(t) and 
width of the kink l{t) which form the basis of their theo- 
retical analysis of our equations of motion and which wc 
repeat below in Eqs.(6) and (8). Where we show their 
equation of motion for P(t) is totally irrelevant to our 
derivations of the SG quasimode and their equation of 
motion for T{t) is independent of the phonon dressing 
whose interaction with x gives rise to the SG quasimode. 
The slope r(i) = 2t: [l{t)]-\ 

We start with the equation of motion for X{t), r{t), 
and x(i) whose solutions are rigorously equivalent to the 
solution of the SG partial differential equation 
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The equation of motion for X{t) is 
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is the bare mass of the kink associated with the X mo- 
tion, and bx = T{a"\x)/Mx- The dot product {f\g) is 
defined as 



The equation of motion for T(t) is 
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Since the center of mass motion does not play any role 
in the existence of the SG quasimode we set X{t) = 
in the equations of motion for x(t) and T{t). We are 
interested in small oscillations of the quasimode so we 
linearize Eq.(5) and Eq.(8) to first order in x E^nd obtain 
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Since we are considering only small oscillations in x '^s 
further linearize Eqs.(lO) and (11) in 6T = F(t) - Fq. 
Finally we obtain 
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due to a potential V{x) is 
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In the remainder of Ref.0] we solved these equations 
of motion analytically and calculated the lifetime of the 
quasimode. 



III. FIRST NUMERICAL INVESTIGATION 

The first numerical search for the SG quasimode in 
Ref.[3| consisted of trying to find the SG quasimode by 
measuring numerically the absorption spectrum of a dis- 
crete SG equation driven by an ac field. What they mea- 
sured was the phonon absorption spectrum of the lin- 
earized discrete SG equation in an ac field which is: 

Xn - iXn+l+Xn-l - 2 X«) + Tq Xn 

= ian+i+an-i-2an)X^ + fit), (14) 

where cr„ is the discrete SG soliton at position n, and Xn 
is the discrete SG phonon at position n. The external ac 
field is f{t) — e exp{iujt) and X{t) is the velocity of the 
center of mass of the SG soliton. We point out again our 
derivations and simulations were for a stationary contin- 
uum SG where X{t) = 0. The spectrum they obtained 
by numerically solving the ac driver discrete SG equation 
is given in their Fig.(l) namely 
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for n^l,2,...,N = /AX, 

(15) 
which is just the spectrum one obtains by solving Eq.(l) 

analytically. For /(t) — ecosujt and X{t) = the spec- 
trum is J2n <5((jJ — (jJ„). If you include X{t) you get ex- 
actly the same spectrum by taking f{t) = cos[(aj/2) t] 
because P^{t) is then proportional to cos[2{uj/2)t], which 
also yields the identical spectrum ^„ 6{uj — uJn) which is 
what they observe and what one obtains by analytically 
solving Eq.(14). Strangely Eq.(14) which is the basis of 
their first numerical investigation, is never mentioned in 
Ref.Q only the numerically observed spectrum Eq.(15) 
is presented. 

What is most important about their first numerical in- 
vestigation is that it could not possibly detect the quasi- 
mode even in principle. In order to observe a quasimode 
in absorption it must first be created and then observed 
during its finite lifetime. Consequently, a quasimode is 
usually observed in emission. The SG quasimode can be 
excited as an initial condition by deforming the slope or 
the width of the kink as an initial condition as we did in 
our derivations and simulations in Ref.[l| and as the au- 
thors of Ref . Q] did in their second numerical investigation 
which we discuss in the next section. The quasimode can 
also be excited by any potential that distorts the slope 
of width of the SG soliton. The force on the slope T{t) 
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For an ac field V — f{t), so the force F vanishes be- 
cause f{t) J dx{da/dT) — 0. Thus an ac field cannot 
possibly excite a phonon mode. Consequently, their first 
numerical investigation in Ref.Q could not possibly de- 
tect the presence of the SG quasimode and thus it has 
no relevance whatsoever to the existence or nonexistence 
of the SG quasimode. It is important to stress that a 
quasimode is different than an eigenmode of a linearized 
Klein-Gordon equation in that an unoccupied eigenmode 
exists even if it is unoccupied. Whereas a quasimode has 
first to be created in order to be observed and it lasts 
only for its lifetime. 



IV. SECOND NUMERICAL INVESTIGATION 

In Ref.jlJ we performed simulations that verified our 
analytic solutions for the SG quasimode. We performed 
three simulations for r(t) and x(0 ^O'' ^ stationary SG 
soliton. We also simulated the Fourier transform of T{t) 
which gives the quasimode frequency and lifetime. The 
simulations agree very closely with analytic results. We 
considered cases where the initial slope was different than 
Fo i.e., Sr{0) ^ and for 5t{0) ^ 0. We considered 
two cases where the length of the system was 1000 units 
and a third case where the length of the system was 200 
units. For the system of length 1000 we followed the time 
development of r(i) and x{t) for times t that were short 
compared with the time a spontaneously emitted phonon 
would travel to the end of the system reflect and inter- 
fere with phonons emitted later. In the third case we 
took a short system L = 200 and followed the time un- 
til the spontaneously emitted phonon first reflected from 
the end of the system. We pointed out that eventually 
the first emitted phonons would reflect form the end and 
interfere with phonons emitted later. Thus simulation of 
the stationary SG quasimode should have a sufhcicntly 
long system that there are no reflections during the time 
of observation or equivalently for a fixed length L the 
time of observation should be less than (L/c) where c is 
the speed of sound. 

In Fig. (2) of Ref.gl the length of the system was L = 
100. They followed the time development of the width 
l{t) for 2500 seconds. The round trip time of a phonon 
emitted by a stationary SG soliton at one end reflect and 
go back to the stationary SG soliton, is 200 seconds. Thus 
during their 2500 second observation time there were 
phonons that made more than twelve trips that would 
interfere with phonons emitted earlier and later. The 
quasimode lifetime is 500 seconds. So phonons could be 
emitted, travel to the end of the system, reflect and be re- 
absorbed by the still excited quasimode. Consequently, 
during their 2500 second simulation time, phonons are 
continuously being emitted, interfering with previously 



emitted phonons reflecting from the ends of the system 
and sometimes being absorbed by the stationary SG soh- 
ton at the end of the system. Consequently, the simula- 
tion of Fig. (2) in Ref.^ should show a very complicated 
interference pattern with multiple time scales but with a 
period of 200 seconds playing a prominent role, which is 
precisely what they observe. If they had taken a much 
longer system or had just simulated for times up to 200 
seconds, they would have verified the existence of the SG 
quasimode. Actually, the first 200 seconds of their simu- 
lation of the width l{t) gives a very good representation of 
the SG quasimode. Thus the incompetent design of their 
simulation in Fig. (2) is the cause of their meaningless, 
complicated interference pattern. 

Once again, in their second numerical solution the au- 
thors of Ref.|2j completely ignore the phonon dressing 
which gives rise to the SG quasimode. In Ref.|l| we cal- 
culated and simulated the phonon dressing. Which shows 
how the dressing decays as the quasimode emits phonons 
during its lifetime while the slope decays from r(0) to 
the constant slope Fq and F(0) decays to zero. Also, as 
in their first numerical simulation, the second numeri- 
cal solution is for an appreciably discrete phonon system 
while our derivations and simulations were for the con- 
tinuum SG equation. Here its interesting to observe that 
the qualitative behavior of the appreciably discrete SG 
system is similar to our analytic calculations and simu- 
lations for the continuum SG. 



In their second numerical investigation of the SG equa- 
tion they started a discrete stationary SG with an initial 
rate of change of the slope f (0) ^ and F(0) 7^ 0. In 
Ref. 1] we considered three such cases except our deriva- 
tions and simulations were for the continuum SG and for 
initial values St{0) of 0.01 and 0.001 and JFq = 0.1. The 
quasimode we derived analytically and verified by sim- 
ulation was for a linear mode. However, they actually 
observed the quasimode in the first 200 seconds of their 
simulation. They however took a system too short for 
the length of time they followed the simulations. Con- 
sequently, they obtained a very complicated phonon in- 
terference pattern due to the multiple phonon interfer- 
ences due to the earlier emitted phonons interfering with 
phonons emitted earlier and later because of the multiple 
reflections of the phonons from the ends of the system. 
In addition, there were multiple absorptions and reemis- 
sions of the phonons with the stationary soliton. If they 
had increased their system from L — 100 to L = 400 and 
followed the simulation for t = 500 seconds instead of 
their t = 2500 seconds they would have obtained essen- 
tially the same diagram we obtained for F(i). 

The authors of Ref.jJI state their theoretical analysis 
of our paper is based on their two cc equations: 
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where 



P{t) = MoloXl-\t), 



V. DISCUSSION 

The authors of Ref. Q] performed two numerical inves- 
tigations for which they state "we show that neither in- 
trinsic internal modes nor "quasimodes" exist in contrast 
to previous reports" referring in particular to our Ref. p] . 
In Sec. Ill we proved that the SG quasimodes that we 
had derived analytically and verifled by simulation could 
not possibly be observed by their flrst numerical inves- 
tigation. The reason is that in order to observe the SG 
quasimode it must first be created and then observed 
during its finite lifetime. We proved in Sec. Ill that an 
ac driver can not create a SG quasimode and thus their 
ac absorption numerical experiment could not possibly 
observe the SG quasimode but could only measure the 
phonon absorption spectra of their discrete SG phonon 
eigenmodes. Our derivations were, for a force free, sta- 
tionary, continuum SG soliton. However, their first nu- 
merical investigation is for an ac driven discrete SG soli- 
ton. They measured numerically the discrete SG spec- 
trum. However, the phonon spectrum of the SG plays 
absolutely no role in our analysis, doesn't appear in any 
of our derivations and is totally irrelevant to our results. 
The SG quasimode comes from the solution of the cou- 
pled continuum equations for the slope of the kink, F(i), 
and the phonon dressing for xit). Consequently, their 
first numerical investigation has no relevance whatsoever 
to the existence or nonexistence of the SG quasimode. 
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Their width variable l{t)/lo is essentially the inverse of 
our variable r(t). Since our X{t) =0 their variable X{t) 
should be identically zero and have no relevance to any of 
our derivations and verifying simulations. They obtained 
the width l{t) in their numerical solution of the discrete 
SG in their Fig. (2) which we discussed in detail in Sec. HI. 
Furthermore, their Eq.(3) for l{t) is incorrect because it 
contains none of the many terms proportional to x(i) that 
appear in the exact equation of motion for T{t) in Eq.(8) 
which are necessary for the existence of the quasimode. 
Consequently, their two equations of motion for P(t) and 
l{t) which they state is the "basis of their theoretical 
analysis" of Ref.Q] have absolutely no relevance to our 
analytic derivation and confirming simulations. 

One of the strangest aspects of Ref. ^] is the complete 
lack of any mention or discussion of the continuum states 
X of the SG equation in the presence of the SG soliton 
that are responsible for the existence of the SG quasi- 
mode. The solution for x derived in Eq.(ll) of Ref.[l| 
constitutes a dynamical dressing of the Sine-Gordon soli- 
ton due to the oscillation of F(t). 

Several times in Ref.[3| the authors compare the SG 
quasimode with the 0^ equation internal mode of the 



which is an exact eigenmode of the hnearized cj)'^ equation 
whose eigenfrequency is in the phonon gap. They state 
that since they can observe the 0"* eigenmode in absorp- 
tion but cannot observe the SG quasimode in absorption, 
this proves the SG quasimode doesn't exist. However, a 
quasimode is not an eigenmode. A quasimode can be ob- 
served in emission but cannot be observed in absorption 
unless it is first created and then observed within its fi- 
nite hfetime. In Ref.[l| we proved that the SG quasimode 
cannot be excited by the ac field used in Ref.0|- 

Finally, in spite of the fact that all our derivations and 
verifying simulations were done for the continuum SG, 
the simulations and analysis by the authors of Ref.y| 
were done for an appreciably discrete SG equation. In 
particular they report finding a discrete mode in the 
phonon gap that was found by Kerekidis and Jones 
Ref. U . It is important to observe that none of the dis- 
crete simulations are in any way relevant to the exact 
analytic calculations and simulations of our continuum 
SG equation. 



VI. CONCLUSION 

In Ref.^] the authors did not make a simple refer- 
ence to or comment upon the exact analytic solutions 
and their verifying simulations of the SG quasimode in 
Ref.[l|. Furthermore, they never mention the continuum 
SG states x that interact with the SG soliton to form 
the SG quasimode. The authors of Ref.i3| performed 
two numerical investigations of the discrete SG equa- 
tion to attempt to prove the SG quasimode doesn't ex- 
ist. The first simulation, an absorption measurement of 
an ac field, which we proved could not possibly observe 
the SG quasimode in principle. The second simulation 
was for an initially deformed discrete SG soliton, which 
they followed in time. In the first 200 seconds they ac- 
tually observed the discrete SG soliton. However, they 
took a system too short for their 2500 seconds observa- 



tion time. Consequently, they observed a complicated 
phonon interference pattern caused by the multiple in- 
terferences between phonons emitted at different times 
which they incorrectly interpreted as the nonexistence of 
the SG quasimode instead of the fact they observed the 
system over twelve times too long for the length of their 
system. 

Recently G. Kalbermann[5J found new important an- 
alytic nonperturbative solutions to the SG equation us- 
ing the Inverse Scattering Transform method. He states 
"his solutions agree very closely with the results of R. 
Boesch and C.R. Willis (our Ref.[l|) in the quasimode 
regime" as shown in his Fig. (4) of Ref.Q- Also, Kalber- 
mann in Ref.j^ points out " a probable source of error 
in the numerical calculations of Quintero et. al.Q-" The 
numerical calculation he refers to in his Ref. |a| is the same 
numerical investigation which we labeled the Second Nu- 
merical Investigation in Sec. IV of the present paper. His 
explanation is essentially identical to the explanation we 
give in Sec. IV of the present paper. 

In their conclusion Quintero et.al. J3| provided an expla- 
nation of our results. They state "the resonance observed 
by Boesch and Willis took place in fact with the lowest 
frequency phonon in the presence of a moving kink and 
not with any internal quasimode." Their argument is, 
since the kink is moving, there is a Doppler shift i.e. 
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C.,^(l + fc2)l/2. 



However, all our analytic calculations and verifying sim- 
ulations are for a stationary kink i.e., X{t) = so the 
kink is not moving and thus the Doppler shift is identi- 
cally zero. Therefore their explanation that our results 
are due to a Doppler shift cannot possibly be correct. 
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